DESY preprint 05-069 



Vortices on the worldsheet of the QCD string 



Harvey B. Meyer 
DESY-Zeuthen, 
Platanen AUee 6 
D- 15738 Zeuthen 

harvey.meyer@desy.de 



Abstract 

We investigate the properties of the QCD string in the Euchdean SU(A^) pure gauge 
theory when the space-time dimensions transverse to it are periodic. From the point of 
view of an effective string theory, the string tension a and the low-energy constants 
of the theory are arbitrary functions of the sizes of the transverse dimensions L± . Since 
the gauge theory is hnearly confining in D = 2, 3 and 4 dimensions, we propose an 
effective string action for the flux-tube energy levels at any choice of L±, given (t{L±) 
and Ck{L±). The Liischer term only depends on the number of massless bosonic degrees 
of freedom and the effective theory can account for its evolution as a function of L± . As 
the size of one transverse dimension is varied, we predict a Kosterlitz-Thouless transition 
of the worldsheet field theory at a{L±)L'^ ~ l/Svr driven by vortices, after which the 
periodic component of the worldsheet displacement vector develops a mass gap and the 
effective central charge drops by one unit. The universal properties of the transition are 
emphasised. 



1 Introduction 



Consider the Euclidean SU(A^) pure gauge theory on a four-dimensional hypertorus of 
size (Li, L2, iv3, L4), with a pair of Polyakov loops running along the direction 4. L3 
and L4 = L are fixed but very large compared to all other length scales in the setup. 
Direction 4 is interpreted as the time direction, so that the potential between two static 
sources can be extracted from the Polyakov loop correlator: 

(P4(0)P4*(^e3)) = e-^(^'^i'^2)L^ (1) 

We are now interested in the system as a function of the size of the two other di- 
mensions, Li and -L2. If one of them is kept large while the other is shrunk to zero, 
we will observe the 'deconfining phase transition', where the centre symmetry along 
the finite dimension breaks spontaneously; decreasing it further, eventually the theory 
dimensionally reduces to a three-dimensional SU(A^) gauge theory. If we now follow the 
same evolution, gradually reducing another dimension, we again go through a deconfin- 
ing phase transition, followed by a dimensional reduction to a two-dimensional theory. 
At this point, the pair of Polyakov considered above loops still obeys an area law; the 
theory is now perturbatively solvable. 

In fact, at any stage of this two-step dimensional reduction, the Polyakov loop cor- 
relator obeys an area law characterised by a string tension cr(Li,L2) varying with the 
size of the transverse dimensions: 

V{R;Li,L2) ~ a{Li,L2)R, La)^^ » 1. (2) 

An effective low-energy theory was recently proposed by Liischer and Weisz ^ to de- 
scribe the corrections to the classical string behaviour V{R) = aR, for the case of infinite 
transverse dimensions in the three- and four-dimensional SU(A^) gauge theories. The 
central idea is to equate the Polyakov loop correlator with the partition function of the 
fluctuating LxR surface, corresponding to the string worldsheet. The leading correction 
to the ground state energy^ comes from Gaussian fluctuations: 

nR) = -R-^^ + o{^). (3) 

where d± is the number of (infinite) transverse space-time dimensions. It is the famous 
Liischer correction jjj , in favour of which a significant amount of numerical evidence has 
been gathered [SI EH . 

Since we have a linearly rising potential at any choice of {Li,L2), we can attempt 
to write an effective string theory for every (Li,L2)- As the latter are varied, we ex- 
pect the effective theory to interpolate between the purely four-dimensional, the three- 
dimensional and the (trivial) two-dimensional case. The unknown parameters of the 
effective string theory now become functions of (Li,L2). 

When the transverse dimensions are large, cr becomes independent of their size. 
As Li shrinks to zero, a oc l/Lf (by asymptotic freedom); when L2 — > as well, 
a is proportional to l/(LiL2). We note that to discuss a smooth evolution of V{R) 
as a function of (^1,^2), we must measure R and V{R) at each (Li,L2) in units of 

^In this paper we shall be interested only in the ground state of the quark-antiquark system, whose 
contribution dominates the correlator at large enough L. 
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a{Li,L2). From the point of view of the four-dimensional theory, the increase of the 
latter as the transverse volume decreases can be interpreted in terms of the overlap 
around the transverse dimensions of the chromoelectric flux lines. Assuming that the 
flux-tube energy density falls off transversely as e""''^^-'- , where mc is the scalar glueball 
mass (something that is hard to establish numerically ^3), the corrections on cr{Li, L2) 
would be 0(e~™'3^™), when L^i^ = min(Li,L2) is large. 

In addition, there are the 'low-energy constants' appearing in the effective action 
describing the energetics of the worldsheet fluctuations beyond the quadratic order. 
These functions of (Li,L2) have to be given as input in the effective theory. The only 
limit in which they are known is the Li,L2 — > limit, where they vanish. 

However, the Liischer correction is completely independent of both a and the c^. 
The question that motivates the present paper is then the following. When the trans- 
verse dimensions Li, L2 are large, the coefficient of the Liischer correction for a bosonic 
string is — 7r/12; when a transverse dimension shrinks to zero, the gauge theory should 
dimensionally reduce to a three-dimensional one, so the Liischer coefficient is then only 
— 7r/24. The 1/R correction disappears altogether when a second dimension is shrunk 
to zero. How and at what temperature do these transitions occur? From the point 
of view of the effective string theory, the coefficient is nothing but the central charge 
of the 2d-conformal theory describing the fluctuations of the worldsheet; it effectively 
counts the number of massless bosonic degrees of freedom living on the worldsheet. We 
thus expect the Liischer correction to have a discontinuity at a phase transition of the 
worldsheet quantum field theory. 

The component along a compact direction (say that of size Li, while L2 stays large) of 
the worldsheet field h{zo, zi) that describes its fluctuations becomes a periodic variable. 
Since the worldsheet field theory is treated as an effective theory with a fixed cutoff, it can 
be described as living on a lattice with a fixed lattice spacing of order a~^^'^. One is then 
led to consider, as effective partition function, the class of statistical mechanics models 
which share the main features of the famous XY- model [3| . In particular the vibrational 
modes of the string are expected to undergo a phase transition of the Kosterlitz-Thouless 
(KT) type |H] related to the proliferation of vortices on the string worldsheet. They 
have the effect of disordering the fluctuations by generating a finite screening length. 
These vortices are interpreted as tunnelling amplitudes between classical configurations 
of the string with different winding numbers. The KT criterion, based on an energy 
vs. entropy argument, yields an estimate of the gauge-theory temperature where the 
transition occurs (see section 4), 



The precise value of this temperature is non-universal; in particular, the string tension 
appearing in the formula is cj(-Li, L2) and differs by a factor of order one from the string 
tension at Li = 00. Since the high-temperature phase of the XY model is dominated 
by vortex configurations and admits a finite mass gap (even in the infinite worldsheet- 
volume limit), this transition naturally explains the disappearance of the universal string 
correction associated with the periodic component of the worldsheet field h. One of the 
aims of this paper is to describe the signatures of such a Kosterlitz-Thouless transition 
in the spectrum of the QCD flux-tube. In this respect, one might worry that T* is 
parametrically of order one in units of the string tension, so that, in the four-dimensional 




(4) 



2 



language, there is no parametric suppression of the string internal degrees of freedom 
at that temperature. But these effects are described in our framework through the 
dependence of a and the on the size of the transverse dimension. 



The rest of this paper is organised as follows. In section 2 we review the properties 
of the 2D XY model relevant for our discussion of the effective worldsheet theory. The 
notation for the latter are set in section 3, and section 4 describes the main results on 
the effect of a finite periodic transverse dimension. A summary and outlook follows. 



2 On the Kosterlitz-Thouless phase transition 

We review some well-known facts about the XY-model and the Kosterlitz-Thouless 
phase transition. For more details, we refer the reader to the abundant literature on the 
subject [Ulllllllin]. 

We take as starting point the 0{2) non-linear cr-model. Once the 2-component field 
is parametrised as (cos 0, sin 0), the continuum action is deceptively simple: 

5 = 11 Sxd^ed^e. (5) 

What makes the dynamics non-trivial is that 9 is an angle variable. Therefore it cannot 
be rescaled, and k is not a redundant variable. The (IR-regularised) two-point function 
in the Gaussian approximation is thus given by — 2mJ^°s(r~)' Such a perturbative 
analysis around the trivial vacuum = of course does not take the periodicity into 
account, and therefore the dynamics is that of an ordinary massless scalar field in that 
approximation. A convenient way to deal with the periodicity of 6 is to let it take values 
in R but identify points differing by 27r. This means that a configuration of can be 
multivalued, as long as e*^ is not. 

The classical equation of motion /S.6 = 0, together with the boundary condition that 
9 goes to a constant at infinity, implies that either 9 is constant or has singularities. If 
the latter are isolated, the most general solution is [TT ^ 
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The interpretation of the solution is that vortices of intensity ±1 are centred at the 
points zf . Higher intensities are obtained when several zf (or several z~) coincide. 
The angle is a multivalued function, but e*^ is well-defined everywhere, except at the 
singular points. 

Although the energy of vortex configurations is divergent both in the infrared and in 
the ultraviolet, in the (lattice) regularised theory their entropy is of the same functional 
form^ log-L/o. So, depending on the temperature, they can play an important role in 
the partition function. The estimate for the critical temperature one obtains from this 
Kosterlitz-Thouless criterion is 

^* = - (7) 
vr 

is the linear size of the system; the energy of a vortex-antivortex pair is roughly obtained by 

replacing L by their separation. 
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These considerations show that to compute the corrections to the free theory, we 
need to formulate it on a lattice and generalise the perturbative approach by taking into 
account Gaussian fluctuations around not only the trivial vacuum, but also the vortex 
configurations: 

S[9 = 0vort + SO] = S^ort [^vort] + f / (d^mf (8) 

Since we have to sum over the different vortex configurations ^vort; the system now 
appears to describe a gas of classical charged particles interacting via a two-dimensional 
Coulomb potential and globally neutral J2 (the neutrality constraint comes from the 
infrared divergence in the 9 two-point function). From one's intuition of a Coulomb 
gas, one infers that these charges are bound into dipoles by a two-dimensional Coulomb 
potential and thus form a dielectric medium at low temperature. At high temperature, 
the charges instead form a plasma which has a finite correlation length due to Debye 
screening. 

The continuum action @ thus has a line of fixed points, > k > k* . To investigate 
the theory at all values of k, one may then choose a lattice action, which has to preserve 
the essential periodicity property of the continuum action and have a local minimum at 
6 = 0. The most general form of the partition function is jl^ 



Z = ^^^^^ Yl £ ^-Jim(x,y))+imix,y)ieixye(y)) ^g^) 

(xy) rn{x,y) 



where the sum is over near neighbours {xy) The function J must be even and have the 
asymptotic behaviour J{m) oc m^, so that the Gaussian action Q is recovered in the 
naive continuum limit. 

The simplest choice is the so-called XY model [H]: 

S,y = K,y ^ [1 - cos{e{x) - e{x + ae^))] . (10) 
x,^l 

This corresponds to J{m) = — log/m(-^xy), where Im is the modified Bessel function. 
At very small K^j it is easy to see by strong coupling expansion J2] that the mass gap 
is given by logl/Xxy When the transition is approached from the high-temperature 
phase, the correlation length in lattice units diverges as 




^ = eoexp , (11) 



where the functional form is universal, while b and -fC*y are not. At very large K^y^ 
the naive continuum limit © is the correct one with k = K^y 

2 ____ 

Another standard partition function, with J{m) = jj^, is the Villain model jllj : 

= n^w^)-^(^+«^M))' (12) 

z{e) = \/2^ Yl expj — ^{e-2Tmf\. (13) 
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Since the asymptotic behaviour of J{m) in the XY-model is m?' /2Ky^y, at large the 
Vihain probabihty distribution becomes identical to the XY probability distribution 
with = i^xy Thus the continuum limit is again ((5} with k = in that regime. In 
general, the vortices decrease the effective k (defined as the k that would give the same 
two-point correlation function in the Gaussian approximation jl4jl. li Rq < a is the 
radius of the vortex-core, the relation is 

— = — + - f:^) . (14) 

Thus the correction is exponentially small at large K^. 

A duality transformation brings the Villain partition function into the following form, 
defined on a dual lattice: 



^ m[x} 

(15) 

where m{x) is a dynamical variable taking integer values. It can be shown j26j to be 
equivalent to the (naively discretised) sine-Gordon model 



(16) 



When the fugacity of vortices is small (this is universally the case on an RG trajectory 
close to the continuum for k > n*), one need keep only the terms m{x) = 0, ±1 in (|15)) 
and the bare lattice parameters are then related by 

13 = 2tt^/K^, ml = Stt^K^ e'^'^^/^ (17) 

What is relevant to our discussion is that the continuum limit of the sine-Gordon model 
at 

P* = ^/8^ (18) 

can be used to describe the Kosterlitz-Thouless transition [T7| [22]. 

The sine-Gordon model is exactly solvable in the continuum by bootstrap meth- 
ods for < < Svr, in which range the potential term is relevant renormalises 
but the coupling (3 does not Remarkably, there is even an exactly solvable (light- 
cone lattice) regularisation of the model [3]. 



3 The effective theory for the QCD string 

We briefly remind the reader of the effective theory for the QCD string 6 . It is in 
principle capable of predicting the splittings (and the degeneracies) between low-lying 
energy levels of the flux-tube, once a few low-energy constants have been determined. 
So it is a low-energy effective theory with a finite UV-cutoff set by the string tension. 
Indeed, at energies of that order, once expects the internal degrees of freedom of the 
string to become excited. 

This theory is specified by a partition function associated with an action living on 
the worldsheet i? x L of the string. The connection with the gauge theory observables is 
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that the Polyakov loop correlation function {Pt{0)* Pt{x)) with |x| = i? is proportional 
to this partition function. The latter may be considered as the partition function of a 
non-renormalisable Euclidean quantum field theory in (1+1) dimensions, on a spatial 
'volume' of linear size R and at temperature T = 1/L; or, alternatively, as a statistical 
mechanics system living on a two-space-dimensional lattice with a lattice spacing of 
order o""^/^ (which has nothing to do with the lattice on which the gauge theory may 
or may not be regulated). 

The field living on the worldsheet {0 < zq < L,0 < zi < R) is a bosonic field h 
with D — 2 components. Since it represents the displacement of the world sheet from 
the classical configuration, it has engineering dimension of length. The Polyakov loops 
are the propagators of static quarks, so we are discussing the open-string case, and this 
is reflected in the effective theory by imposing Dirichlet boundary conditions on h at 
zi = and zi = R. The action of the effective theory is 

S = aRL + fiL + So + Si + S^^^ + S^^^ + ... (19) 
where the Gaussian action 

So = I y <fz d^h ■ d^h (20) 

gives the leading contribution to the string spectrum at large R. Each higher-order 
term is multiplied by an unknown dimensionless constant , times a; some of them 
vanish/are constrained by the open-closed string duality [S]. The lower index on the 
action terms above indicate how many more derivatives of h they contain than Sq. The 
parameter a controls the amplitude of the fluctuations around the classical solution: in 
this form, the whole action of the effective theory is multiplied by it. Of course, a is 
of the same order as a itself, and the fact that we are to treat the latter as large (it 
is the UV cutoff) shows that we are in fact doing a semi-classical calculation. But if 
we introduce the canonically normalised field (p = \pdh^ a is absorbed into the low- 
energy constants, = j ia)^!"^ . The value of oja is not observable here: if for 
instance the low- lying energy levels are determined from lattice simulations, one will 
simply determine the dimension- full constants in units of a. 

It should be noted that for any relativistic bosonic string, one has in fact a = a: 
the amplitude of Gaussian fluctuations are dictated by the same string tension that 
determines the energy of the classical configuration. It appears however that one cannot 
test this property for the fiux-tube by measuring the low-energy spectrum in infinite 
transverse dimensions. From now on we drop the tilde on a, but it is good to keep the 
distinction in mind. 

4 The effect of a periodic transverse dimension 

We now consider the effective theory describing the Polyakov loop correlator, for sim- 
plicity, in D = 3 dimensions - we will later return to the D = 4 situation. The transverse 
dimension (the one in which the unique component of h can fluctuate) is periodic, with 
period Li = L±: 

(P3(0)P3*(«e2))^e-^(«'^^)^, L3^L»i?. (21) 



6 



How is one to generalise the effective string theory discussed in the previous section 
to include the effect of the periodicity of the transverse dimension? 

The general requirements on the effective action (described in detail in ^) that lead 
to the form (|2()|) must be kept, but h is now a periodic variable: h+L± must be identified 
with h. It is then convenient to rescale h according to 

h = (22) 

Because of this periodicity, a convenient way to write the most general partition function 
is to formulate it on a lattice with lattice spacing of order 1 / \/cf{L_[_). Working on a 
lattice is equivalent to working with a continuum effective action with a fixed momentum 
cutoff |23| . Our Ansatz for the effective partition function is 

Z{L^,L,R) = c''^L,)RT+nT J -Q dejx) -Q ^ ^-j{L^M^,y))+im{^,ym^)-e{y)) 

^ (xj/) m{x,y)&Z 

J{LjL,m) ~ mV(2K), m ^ oo. (23) 



where (xy) is a pair of nearest neighbours and 

(2vr) 



K^-^i^. (24) 



The naive continuum limit is then given by 



So = yJ ^A*^ • (25) 

Interactions over a few lattice spacings are possible but should not affect the universal 
properties of the system, as long as the naive continuum limit is the one given in (|25|) . 
The simple application of the Kosterlitz-Thouless analysis leads to the following scenario 
function of L± (assuming R,L very large): 

• for large L± (low-temperature phase of the XY model), the Gaussian approxima- 
tion is good, and the periodicity of 6 plays no role for the long-distance properties 
of the theory. In particular, the mass gap vanishes and the universal part of the 
Casimir energy of the system remains exactly 

• as L± decreases, the vortices can affect the non-universal contributions to the 
energy levels. 

• when L_|_ reaches L*j_, vortices on the worldsheet become relevant and the Kosterlitz- 
Thouless transition takes place. The KT estimate of L*, is given by 



for L± < L*j_, the theory has a finite mass gap m. The Casimir energy if of order 
g-2mL^ and the fluctuation along this dimension is effectively 'switched off'. 

When L± <C L*j_, the strong coupling expansion (see e.g. ,12^) tells us that to 
leading order in k, m = — log k. 
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It can be shown |14j that the KT transition takes place at somewhat larger L±_ than 
the KT estimate because of the disordering effect of vortices. K* can be extracted 
from lattice simulations: first the function ct(L_|_) must be computed, and then K* can 
be obtained from the determination of L\. It should be noted that it is a non-universal 
quantity; for the square-lattice XY-model, it was found in numerical simulations [TSj that 
AC* = 1.1199(1), relatively far from the Kosterlitz-Thouless estimate of 2/tt ~ 0.6366. 

Physical interpretation of vortex configurations It is useful to complement the 
machinery of two-dimensional statistical mechanics with some more intuitive consider- 
ations. The action (|23|) is a saddle-point expansion around the classical string solution 
which connects a quark to an antiquark separated by distance R. Clearly, there are also 
classical solutions with a net winding number around the compact transverse dimension. 
They have energy 

Sd(n) = a(L^)/R^T^=.i?(n = 0) + ^^:^^i^, fl»L^. (27) 

The n = 1 classical solution becomes as light as the first Gaussian excitation of the 
n = solution when cr(L_|_)L^ = 27r. The different classical configurations can only 
be connected if the string goes through more energetic configurations. Therefore such 
transitions are classically forbidden tunnelling processes. Quantum mechanically, one 
expects the tunnelling to become a frequent fluctuation when the zero-point energy 
of the n = vacuum is of same order as the gap E{1) — E{0). In this context, the 
worldsheet direction of size R is interpreted as the space-direction, and that of size L 
as the Euclidean time direction in a path integral treatment of the quantum-mechanical 
string. A vortex being a point-like object on the worldsheet, it corresponds to a process 
of the string describing a transition from a state at t = — oo to another state at t = +oo. 
As long as we are considering the asymptotic limit L — > oo, only vortex-antivortex 
configurations have a finite free energy in the 2d QFT. By writing the solution @ 
with a single vv pair explicitly, it is easy to see that a worldsheet containing a vv pair 
describes the 'life' of a string that had winding number n = at t = — oo, goes through 
a state with winding number ±1, and returns to n = classical state (see Fig. 

One may find the discussion analogous to the discussion of instantons in the 4D 
infinite volume gauge theory. The (anti)vortices correspond to (anti)instantons, the 
classical configurations of the string correspond to the configurations of the gauge field 
at fixed t where the mapping from the sphere S'^ at spatial infinity to SU(2) has a 
deflnite winding number. A difference is that, due to infrared divergences, in the 2d 
system only configurations satisfying J2m "^(™m ~ ^m) = are allowed, where rim (^m) 
is the number of (anti)vortices of strength m in the configuration. By contrast, there 
are configurations of the 4D gauge theory with a net topological charge which have a 
finite Euclidean action. 

EfTect of the vortices on the string spectrum at L± > L*j_ For L± > L*j_, we are 
in the low-temperature phase of the XY model, and vortices correspond to an irrelevant 
operator. Therefore they do not affect the universal 1/R correction, which is determined 
by the canonical kinetic term, the only renormalisable one compatible of the symmetry 
requirements on the effective action. 
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Figure 1: Physical interpretation of the vortices: the left-hand side, where 9 is repre- 
sented as a unit vector (cos 6, sin 9), illustrates that the appearance of a vortex at time t 
'followed' by an antivortex at time t induces one unit of clockwise winding of the variable 
9{ts,x) as X varies from to i? on a timeslice at time tg located between t and t. Right, 
the corresponding interpretation for the string on the covering space of the torus: the 
string has non-trivial winding number along the transverse dimension in the interval 
[t, t\ and trivial winding number elsewhere. 
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More generally, Ec{R)R and Ei(R)R, the Casimir energy and the energy gap to 
the first string excitation, are low-energy physical observables in the fixed-cutoff world- 
sheet quantum field theory; in fact, the latter quantity was used as a definition of the 
renormalised coupling in the XY-model in (2^. According to Symanzik's theory on 
the approach to the continuum limit of physical observables, irrelevant operators give 
a contribution which is suppressed by the lattice spacing to the power of their scaling 
dimension minus the number of space-time dimensions j23|. In particular the higher- 
derivative interactions on the worldsheet give contributions to the energy levels which 
are suppressed with respect to the 1/R term by 1/ \/(t{L±) to the power of the num- 
ber of derivatives of the irrelevant operator, minus two. Now as long as the theory 
is scale invariant, dimensional analysis tells us that the power of 1/R that suppresses 
the contribution of an irrelevant operator to E{R) ■ R must be this same power. For 
instance, (d^cp)'^ {8^4')'^ contributes as l/{a{L±)R'^) to E ■ R, if E is a string energy 
level. Thus the parametric form of the energy levels computed in i6_ could be guessed 
on general grounds, while relating the coefficients to the action and computing the 
string degeneracies requires an explicit calculation. Similarly the vortex-operator gives 
a suppressed contribution. Its scaling dimension, as computed in the Gaussian theory 
S= ^j(fx{d^ef, is [El 

= VTK. (28) 

Now, since the free energy of a single vortex is divergent in the infrared, the single- 
vortex contribution to the partition function vanishes in the thermodynamic limit of the 
2d system jl5| . The first contribution comes from a vortex — anti- vortex pair, so that 
the vortex contribution to the string ground state energy is suppressed by \/R to the 
power 2x^ — 2: 

AE{L^,R) = ^{a{L^)R')'--^ (29) 

The peculiarity is that the power varies continuously with the transverse size. We note 
that it can also be obtained by using the connection with the sine-Gordon model and 
the fact that the scaling dimension of the cos (34>{x) operator is ^ . 

We now discuss the relation between k and L± . The lattice theory, at lattice spacing 
a = 1/ y/a{L±), is initially parametrised by the coupling K = a{L±)L'j_/ Att"^ . Under 
the renormalisation group the coupling K evolves as a function of fia, where ;U is a 
renormalisation scale, in order to keep the long-distance physics constant. Under coarse- 
graining, the theory flows into a Gaussian fixed point, with a Gaussian coupling k = 
lim^a^oo K{^a), which is a non-universal function of the initial K. All we can expect is 
that K is monotonously increasing in K. 

To determine the RG flow of -R' in a given model, we have to choose a physical 
quantity that we keep flxed as the lattice spacing is varied. One possibility is to choose 
the "spin-spin" correlation function, (e*^^'-''-*"^^^'-*^), which in the Gaussian theory behaves 
as C{a/\x\)^^'^'^'^ at large distances. In the lattice theory at lattice spacing a, the long 
distance behaviour of this correlator will also decay algebraically, with a certain exponent 
b. Then the renormalised K{na) can be set to l/27r6. In the Villain (also called 'periodic 
Gaussian') model, this computation can be carried out explicitly (see section 2, Eg. ITU) . 
At large K = K^, the shift from K to K{^a) is exponentially small in K. In that speciflc 
model, the vortices are the only origin of renormalisation. The exponential suppression 
of their effect is a universal fact, because their fugacity is of order y = e~'^^ , where c is 
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however a non-universal quantity representing the core energy of a vortex. In general, 
there is in addition the more usual renormalisation of K due to operators containing a 
higher number of derivatives. These effects are of order one. So although the scaling 
dimension of the vortex operator is a non-universal quantity, at large K, we have 

ti = K +0(1). (30) 



Properties of the worldsheet phase transition We now proceed to discuss the 
universal signatures of the KT transition, from the point of view of string observables. 

A first point to note is that the critical value of k* = 2/7r is universal (while that of 
K* is not). Therefore the vortex-antivortex pair contribution to the string energy levels 
reaches 

/^E (xl/a{L\)R^ (31) 

at the KT transition. At the KT point, the fluctuation-field correlations on the world- 
sheet read, for ||x|| » l/^o"(L_|_), 

(MO)Mx)) = -(^)'log|N|. (32) 

However the vortex operator becomes relevant at that point, and (anti)vortices start to 
proliferate. 

If one is approaching L\_ from below, that is, from the 'high temperature' phase 
of the XY model, another universal property of the KT phase transition is the rapidly 
vanishing mass gap 

m~moexp . (33) 

m-o, K* and h are non-universal quantities; but given the function cr(L_|_), this translates 
through Eq. [^into a robust prediction for the functional form of m as a function of L±_ 
close to the phase transition. {a{L±^ is defined by the asymptotic static force; we are 
assuming here that (t(L_l) also governs the quantum fluctuations (see section 3).) So 
the KT nature of the transition can in principle be tested, but we note that in practice 
the range where the scaling H33|) holds is typically very small p[5j . 

When a mass gap m appears, the Casimir energy must decrease exponentially in mL 
when the product is large. For a theory of free, neutral massive bosons with Dirichlet 
boundary conditions, one can indeed show (see Appendix A) that 

Ec{L) = + 0(6"^"^^) (34) 

4 V ttL 

For the current theory, which contains vortices and anti- vortices, one expects twice this 
expression; and to obtain the corresponding expression for periodic boundary conditions, 
one can apply the rule PO] 



-E'c(-^^)periodic = ^ £'c ( — ) . (35) 

^ / Dirichlet 



The Casimir energy is thus smaller than any power of 1/-L, and the effective central 
charge vanishes. So the phase transition of the effective theory driven by vortices explains 
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the disappearance of the Liischer correction as a transverse dimension is reduced in size. 
Although the vortices are not free, in the sine-Gordon model they correspond to solitons, 
which are the only one-particle states for > 47r, where ~ a{L±)L'^. Their S-matrix 
is known exactly , and the Casimir energy can be obtained from an exact non- linear 
integral equation [3]: in a periodic box, it reads 

Ec{L) = -—K,{mL) ■ ' ^^"^ 



+ 



TT 



mL 



+ 



1 



(36) 



The leading term is identical with the free, complex, scalar field with the same periodic 
boundary conditions. Thus this analytic form is universal, for a given mass gap m. 



The spread of the cross-over for a finite-length string Strictly speaking, the 
worldsheet statistical mechanics system cannot have a phase transition for a finite-length 
string, instead the latter is generically smoothed out into a cross-over. Decreasing L± 
a little below L*j_ leads to a suppression of the 1/R term, which becomes exponential in 
moR within 

— = (37) 
K* log^moR ' 

From then on, this transverse dimension is effectively "switched off" from the string 
dynamics. The argument above shows that the width in K (and presumably also in L±) 
of the smoothed-out phase transition is only logarithmically suppressed in the string 
length, a distinctive signature of the KT nature of the transition. 

Since mo is expected to be of order \/a{L±), what we just learnt has the practical 
consequence that R must be exponentially large in the correlation length of the under- 
lying gauge theory, for this cross-over to be observed at all. This might well explain why 
these effects were not observed in the numerical studies of spatial torelon masses 

in 2+1 dimensional gauge theories. 



The transition region for a finite-length string No matter how long the string 
is, for L_L sufficiently close to L*j_ the mass gap becomes much smaller than 1/R. To 
get an idea of how the Casimir energy of the string then looks like, we compute in the 
appendix the leading corrections in the small mR regime to the conformal field theory 
value — 7r/24i? in the neutral free scalar field case: 

Ec = --^ + 0{m) + o(m'^Rlog{mR)) . {a-^/^{L^) <. R ^ m-^) (38) 
24x1 ^ ^ 

More precisely, in terms of the static force, we have 

2 

\F{L^,R)\=a{L^) + ^ + '^logmR + ... {a-^'\L^) <^ R <^ m'^) (39) 

It is effectively as if the string tension had a logarithmic dependence on R in this regime. 

In the case of periodic boundary conditions, the leading term in Ec would be — 
This is also the outcome of an exact calculation in the sine-Gordon theory jKj. Inter- 
estingly, in this regime the solitons and anti-solitons do not lead to a doubling of the 
effective central charge; on short distance scales, fluctuations are of the "spin wave type" 
dictated by the kinetic term of a real scalar field. The Casimir energy of the sine-Gordon 
model thus interpolates between the Casimir energy of a free complex scalar field in the 
infrared and a free real scalar field in the ultraviolet. 
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4.1 Additional transverse dimensions & applications 

So far we discussed the dynamics of a single transverse dimension, in other words we were 
deaUng with a three-dimensional gauge theory. Now we return to the four-dimensional 
case, where the two components of the worldsheet fluctuation field h are periodic, the 
case of relevance to standard numerical simulations. 

The simple but key point to note is that the two components are decoupled in the 
quadratic, renormalisable action (|2nj) responsible for the Liischer term. Therefore the 
universal contributions to the string energy levels just add up. Also the existence of 
topological configurations in one component is unaffected by the presence of the other 
component. In general, the two XY-type models for the two components are coupled 
by 'irrelevant' terms. This implies in particular that the energy of a vortex-antivortex 
'dipole' of size L^^ receives contributions of order l/(cr(L_L)L^^). This however does not 
affect the Kosterlitz-Thouless criterion for the occurrence of a phase transition driven 
by the vortices. We therefore do not expect the KT nature of the worldsheet phase 
transition to be affected by an additional fluctuation component. Thus, if starting from 
the 4D theory, we successively reduce two dimensions, the latter get in turn 'switched 
off' (from the point of view of contributing to the central charge) at two different KT 
phase transitions. If we reduce them simultaneously, there presumably are two nearby 
transitions; this situation however deserves a more detailed investigation. 

A convenient way to locate the worldsheet KT transition in Monte-Carlo simulations 
of the lattice gauge theory by working at a relatively small Li = L± < L\ is to measure 
the static quark potential off-axis, say at an angle a: {P4{0) P^{R cos ae^ + Rsinaei)). 
The effective periodicity of the worldsheet fluctuation field is then given by the substi- 
tution L_|_ — > L±{a) = L±/ cosa. By measuring the static potential Va{R), one can 
in principle locate the a* where the Liischer term drops by a unit of central charge. 
One can then check in an independent simulation that the transition is also seen when 
setting L_L = L±{a*) and measuring the static potential on-axis. This would test that 
the transition is really driven by the periodicity of the fluctuation field h{z). 

5 Summary and outlook 

We have studied an effective string theory for the flux-tube in SU(A'") gauge theories in 
three and four compact dimensions. Let us summarise the outcome of our investigation 
of the Polyakov loop correlator when one transverse dimension is periodic and of size 
L±. The periodic component of the fluctuation field is described by an action of the 
XY-model type. It follows that 

• there is a phase transition of the worldsheet field theory at L*j_ = 0(o"^^/^) where 
the periodic component acquires a mass gap and hence the central charge drops 
by one unit; 

• it is a Kosterlitz-Thouless phase transition: 

• approaching L*j_ from below, the mass gap m{L±) goes to zero in a universal way 
in terms of K = a{L±)L'^/A7r'^ as given by Eg. 1331 Once <j{L±) is independently 
determined, we have a prediction for the functional form of ■m{L±). 
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• in the large L± phase, there are (strongly suppressed) corrections to the string en- 
ergy levels associated with the periodicity; they come in powers of which vary 
continuously and increase monotonously with L±; these corrections are exactly of 
order l/R^ at L± = L]_. 

We conclude by mentioning some natural extensions of the investigations carried out 
in this paper. 



1. It is also possible to study the QCD string when one transverse dimension is finite 
with Dirichlet boundary condition for the gauge field, as in the Schrodinger functional 
(see for instance plj ) . The first question to address is what this boundary condition 
translates into for the worldsheet degrees of freedom; presumably it introduces a cutoff 
in the amplitude of the fluctuations. Unlike periodic boundary conditions, this cutoff 
will affect even the "spin waves" of the worldsheet. It is known from the discussion of 
the roughening phenomenon in lattice QCD that the bosonic string wave functional is 
Gaussian in the transverse coordinates and that its width grows with the string length 
i? as 122] 

5\R) = — log ^. (40) 

One would thus expect that until R is 0{e^^^), the effect on the string spectrum 
is suppressed by e~^^^^ At fixed L±, this offers a chance to verify the quantum- 
mechanical broadening of the string. 

2. In string theory the duality transformation L± <-> plays an important role. It 
is clear that a self-dual point at L±_ ~ O.lSfm does not exist in the SU(A^) theory at 
finite N. From the point of view taken in this paper, it is because the string tension 
and the effective string action varies with L±, unlike a fundamental string theory. 



3. The fact that the QCD flux-tube, when probed at energies well below the scale set 
by its tension, behaves like a bosonic string living in the same number of dimensions 
as the gauge theory, is highly remarkable. One might call this the "weak gauge-string 
equivalence" Besides this line of investigation, there is the long-standing speculation 
that the gauge theory itself has an equivalent formulation in string theory (the "strong 
gauge-string equivalence" _27J. Since the gauge theory is renormalisable, it is a stand- 
alone theory which describes phenomena on all scales. Therefore the equivalent string 
formulation must also make sense at all energy scales, which means that it must be 
a fundamental string theory. All known fundamental string theories live in a higher 
number of dimensions. 

It is natural to ask what the relation is between the effective four-dimensional string, 
whose properties can be established through the weak gauge-string correspondence, and 
the hypothetical, presumably higher-dimensional string theory which behaves like the 
gauge theory on all scales. One would naively expect the former to arise as the low- 
energy effective string obtained by dimensionally reducing the degrees of freedom of 
the latter, fundamental string theory. If such a relation exists, the worldsheet phase 
transition of the effective string described in this paper must also affect the fundamental 
string dynamics, and hence should be reflected in other physical observables of the 
gauge theory. The gauge theory of course has the 'deconfining' phase transition at 
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L_L = 1/Tc = 0((T~"^/^), and it is clearly interesting to see whether L\_ = l/Tc- This 
can be achieved straightforwardly by lattice simulations. Since dimensional reduction 
in the gauge theory does not work all the way down to Tc, there is no a priori reason 
to expect L*j_ = l/Tc- such a coincidence would provide evidence for a relation of the 
effective QCD string to a fundamental string. 

The low-energy constants determined in the weak gauge-string correspondence also 
provide valuable information on the 'mother' string theory, and the geometry of the 
space which leads to its dimensional reduction. Their dependence on the number of 
colours is of particular interest [7j, since the equivalence to a string theory is expected 
to be simpler at large N. 



4. We note that the Kosterlitz-Thouless worldsheet transition was discussed in the con- 
text of the bosonic string, as a fundamental theory in 26 dimensions, by Sathiapalan 
The fact that the 0(2) nonlinear cr-model has no fix points for k < n* implies that it 
cannot be the worldsheet action of a fundamental string theory. The author concludes 
that in this regime, the string can no longer be considered to be in a flat background, 
because the vortex condensate contributes to the stress-energy tensor and hence to the 
space-time curvature. In the case of the QCD string, we have a fixed cutoff, and we have 
shown that the net effect is a change of the central charge, as measured through the 
Casimir energy of the string: one component of the worldsheet field effectively acquires 
a mass. Thus an attractive scenario for a fundamental string is the following. When 
probed at low energies E <^ T* , the Casimir energy appears to be that of a string in 
lower dimensions, due to the curvature of the target space-time which gives an effective 
mass to some components of the worldsheet field. But at higher energies, the curvature 
is irrelevant and the "hidden" components can be excited: the original conformal invari- 
ance is restored in this regime. In this way, it is plausible that there exists a fundamental 
string theory which at low-energies describes the QCD string. 
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A The Casimir effect for a massive scalar field 

The Casimir energy for a free massive scalar field in d space dimensions with Dirichlet 
boundary condition in one dimension of size R reads |2nj 



1 / 

Ec(m,R) = ^ / dt f^-^ log (1 



This is the energy per unit of transverse volume. A term proportional to R and a term 
proportional to m have already been discarded at this stage [20]. By expanding the 
logarithm log(2;) around z = 1 and integrating term by term, one obtains Ec as a series 
of modified Bessel functions: 

-2 rmR\ i^) ^ Ka+,{2mRn) 
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Given the asymptotic large-argument behaviour of the Bessel functions, 

this expression of Ec immediately shows that the latter vanishes exponentially fast in 
mR once this product is large: 

Ec{mR -^oo,R)^ (if )'''^"'"''- (43) 

On the other hand, when 

e = mR < 1, 

one naturally expects to recover the Casimir energy of a massless scalar field [20] 

Ec{m = 0,R) = ^ r (^"j C{d + 1). (44) 



This is easily checked to be the case, either by using the small-argument expansion of 
the Bessel functions, 

or by contour integration techniques of (|4H) (with mR = 0). The leading corrections to 
this was given in PUI for d > 3: 

Ec{m,R) = Ec{m = 0,R) + —l^l^T(^^y{d-l)imR)^ + ... (46) 

For d = 2, the leading correction is easily computed by the change of variables = s^— 
in (|ni): 

2 

Ec{m, R) = Ec{m = 0,R)-^ [log {2mR) - 1] + 0{mRf, d = 2. (47) 

We are interested in the Casimir energy of a string, so that d = 1 and, using Eq. 
and C(2) = we have the well-known result |^ 

Ec{m = 0,R) = -^. (48) 

We now want to compute the leading corrections in e to this result, which corresponds 
to the conformal theory. This is a less straightforward calculation, because one cannot 
simply differentiate (with respect to e) the integrand in the first form H41|) . or the terms 
of the series in the second form (|42|) . because in both cases the result would diverge. 
We start off from H42|) . The key point to note is that as e — > 0, the series samples the 
Bessel function in very small steps, so that in can be approximated by an integral: In 
general, for a smooth function g{x), one would have 



gi^n) = / g{x)dx + 0{e^) 

n>l •^^/2 



(49) 
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This is however not directly applicable, because the Ki{x) function is singular at the 
origin. We therefore differentiate Eq with respect to e at fixed R: 



-2.r'-^4^ = E^^ + 2ei.((2en) (50) 
= -2eJ2Koi2en) (51) 

n>l 

where K[ is the derivative of Ki with respect to its argument. We have used the identity 
K[ + Ki/z = —Kq. We now use the formula, valid for an analytic function, 

2e E men) = H /(x)dx - 2e ^ E QtW'"" ' ^''^ 

n>l n>lfc>l + 

Now, with / = Ko, 

/ Koix)dx= Ko{x)dx- Ko(a;)dx = - + e[loge-l-log2 + 7]+0(e^loge). 

Jt Jo Jo 2 

(53) 

The remaining sums in (|52j) can be evaluated in leading order in e by using the leading 
small-argument behaviour of i^o(e) = — loge + 0(1): 

J: K^^'\2ne)e^' = ^'^"^^^^ + 0{eHoge) (54) 

ra>l 

Gathering the pieces and performing a few elementary integrals, we get 

dEc{e',R) 



Ec{m,R) = Ec{0,R) + J^'-^^^de' 



vr in m?R 
+ -r + 



2AR 4 47r 



log — 



^^-4^fc(2A: + l) 

+0(m^i?3 log mfl) (55) 

It is clear from the results above that the introduction of a small mass for the scalar 
field has a stronger effect in lower space dimensions. This is natural since infrared 
effects tend to be stronger in lower-dimensional field theories. The term proportional to 
m?R\ogmR is the first term of interest, in the sense that it is the leading observable 
effect in a system with fixed m and varying R. It expresses the fact that the bulk free 
energy of the system has a logarithmic dependence on its size, in the regime mR < 1. 

The cross-over between the small mR and large mR regimes is smooth. The area 
—R /g°° deEc{e,R) is easily computed and amounts to C(3)/16. 
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